Evenly distributed unitaries: on the structure of unitary designs 
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We clarify the mathematical structure underlying unitary t-designs. These are sets of unitary matrices, evenly 
distributed in the sense that the average of any t-th order polynomial over the design equals the average over 
the entire unitary group. We present a simple necessary and sufficient criterion for deciding if a set of matrices 
constitutes a design. Lower bounds for the number of elements of 2-designs are derived. We show how to 
turn mutually unbiased bases into approximate 2-designs whose cardinality is optimal in leading order. Designs 
of higher order are discussed and an example of a unitary 5-design is presented. We comment on the relation 
between unitary and spherical designs and outline methods for finding designs numerically or by searching 
character tables of finite groups. Further, we sketch connections to problems in linear optics and questions 
regarding typical entanglement. 



I. INTRODUCTION 

Before introducing the notion of a unitary design it is 
worthwhile to look at the analogue structure on spheres in 
IR™. Imagine one is interested in the average value of a real 
function / defined on an n-dimensional real sphere S n . That 
value might be hard to compute in general so it could be sensi- 
ble to estimate it by averaging over a finite set of unit vectors 
V = {\i/ji), ... , a')}- Of course, for any such finite set, 
there are functions whose true average value deviates arbitrar- 
ily much from the one approximated by summing over T>; but 
the more points the test set includes and the more "even" these 
vectors are distributed, the more "exotic" such functions have 
to be. The following notion aims to quantitatively capture the 
quality of a set of points for these purposes: a finite subset 
T> of S n is called a spherical t-design if the average of ev- 
ery t-th order polynomial p over S n equals p's average taken 
over T>. A large body of literature has been devoted to the 
construction and exploration of designs. Many of the relevant 
references can be found in the accessible article Ref. yfl. 

One can adapt the definition of spherical designs to com- 
plex vector spaces (simply by substituting the real sphere by 
the set of complex unit vectors) with obvious applications 
in quantum mechanics. In the context of quantum informa- 
tion theory, 2-designs appeared in Refs. Jl, H, H, 0, H, IH, 
to name a few. Two prominent examples of complex spher- 
ical 2-designs are here known by the names of mutually unbi- 
ased bases 1 3, 8, 10] and symmetric informationally complete 
POVMs d respectively. 

Quite recently, Dankert et al. introduced the notion of a uni- 
tary t-design by replacing the real sphere S n by the set of uni- 
tary matrices 11(d) in the definition of spherical designs [1,2]. 
Averages are here to be taken with respect to the Haar mea- 
sure. In the sense made precise above, the theory of unitary 
designs thus aims to identify finite nets of unitaries, which 
cover the entire group as tightly as possible. 

Such nets are interesting for various reasons. Abstractly, 
unitary designs can serve as testbeds for examining conjec- 



tures concerning the unitary group: their even distribution 
here means that they "cover many complementary aspects" 
of U(d). Also, spherical designs naturally appear in optimal 
solutions to several physical problems, ranging from quantum 
state tomography and optimal estimation using finite ensem- 
bles to quantum key distribution IdI la, 171 1 1111 - it is sensible to 
assume that similar applications for their unitary counterpart 
can be found. More concretely, unitary designs have been ap- 
plied to quantum process estimation and fidelity estimation of 
channels using random states IH Bl; quantum cryptography 
lfl2ll and data hiding protocols {Ofl . Naturally, designs can be 
used to estimate Haar averages using classical computers. For 
the case of averages of polynomial functions, the results are 
guaranteed to be correct (there are, however, other methods 
for tackling this specific problem; see Appendix IVIII Al and 
Ref. lll4ll ). For non-polynomial functions one obtains at least 
an educated guess. Going beyond finite-dimensional quantum 
systems, Haar averages over U (d) appear in the context of 
energy-preserving transformations of d bosonic modes. Such 
transformations are notably relevant as passive linear optical 
transformations of states of light modes IT151 [l6ll . Lastly, we 
believe the problem to be of inherent geometrical interest. 

As most of the present work is concerned with unitary 2- 
designs, we now state the precise definition for this special 
case (see, however, Section [V At : 

Definition 1 (Unitary design HH). A set V = {U k }k=i,...,K 
of unitary matrices on H. — <C d is a unitary 2-design it 
fulfills the equivalent conditions: 

1. (Averages) Let p be a polynomial in 2d 2 variables. We 
can conceive p as a function on U(d) by evaluating it on 
the matrix elements and their complex conjugates of a 
given matrix: p(U) := p{U l j,U l j). One now demands 
that for any p which is homogeneous of degree two in 
each variable, the relation 



~ E P<Pk) = I : 
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be fulfilled. 
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2. (Twirling of states) For all p e B(H <g> H) 

^ ( U k® U k ) p (U k (8 U k )l (2) 
u k ev 

= [ (U ®U)p{U ®U)UU. (3) 

Ju(d) 

3. (Twirling of channels) For any quantum channel A 

i ^ utA(U k pUt)U k (4) 

= / U^K{UpU^)UdU. 

Ju(d) 

The problem has a long history, which is formulated mostly 
in the second of the three equivalent guises listed above. The 
"twirling" operation originates from invariant theory (where 
it is sometimes called "transfer homomorphism") and has, to 
our knowledge, first been introduced to quantum information 
theory in Ref. f\Ja . giving rise to the concept of a "Werner 
state". Later, it was noted that in d = 2 (i.e., for single qubits), 
it suffices to average over a finite set of unitaries iflSfl . A 
construction for general dimensions - employing non-evenly 
weighted unitaries - appeared in Re f. 11911. DiVincenzo et al. 
lfl3ll realized that the Clifford group ll20ll23[l for qubit systems 
exhibits the property given in Eq. (O; a fact which was later 
gen eralized to systems of prime-power dimensions by Chau 
01211 . Similar ideas appeared in Ref. 12111 . A first concise 
treatment was given in a master thesis by Dankert flj] (where 
the term of a unitary t-design has been coined) and in a later 
paper by Dankert et al. fH]. In these publications, the equiva- 
lence of the criteria in Definition[T]has been made explicit and 
the question of how to efficiently implement the unitaries of 
certain designs was addressed. 

Despite the large amount of interest paid to the problem, 
the following natural questions have been left open and will 
partly be answered in this paper: 

1. In which dimensions do unitary 2-designs exist and 
when can they be explicitly constructed? While we do 
not have a general answer to this question, a host of 
examples is provided in Sections [ill] and IIV CI [Note 
added in revised version: After this article had been 
submitted, A. Scott made us aware of Ref. j45tl . This 
extremely general paper proves - among other things 
- the existence of unitary designs for every t and d (it 
does not provide an explicitly way for constructing the 
designs). Thus, the question posed above can partly be 
answered affirmatively.] 

2. What is the minimal number of elements needed for a 
2-design ? See Section III CI for a lower bound, which 
we conjecture to be tight in leading order. 

3. Is there an easy criterion to decide whether a given set 
of matrices constitutes a design? This question is an- 
swered affirmatively in Section III B I We transfer the 
concept of a frame potential yl l3l|] from spherical to 




(a) (b) 

FIG. 1: Visualization of the 12-element Clifford 2-design described 
in Section lTV C| As up to phases SU (2) ~ 5*0(3), every qubit uni- 
tary corresponds to a three-dimensional rotation. The group SO(3), 
in turn, can be pictured as a ball with radius n, where antipodes on the 
boundary are identified. This is done by associating to every rotation 
by an angle <t> 6 [0, 7r] about the unit- vector h the point <ph G K 3 . 
Figure (a) shows the four Pauli matrices 1, a x , a v , a z in this rep- 
resentation. The non-trivial Pauli operations lie on the boundary of 
the ball and hence appear twice: a x , e.g., at ±(7r,0,0) T . Adding 
eight further Clifford operations, which correspond to the vertices 
27r/\/27(±l, ±1, ±1) T of a cube, we arrive at the 2-design pictured 
in Figure (b). 



unitary designs. The frame potential is a simple poly- 
nomial expression in the matrix elements, which is min- 
imized exactly for designs. This criterion even allows 
for numerical searches in spaces of small dimensions. 

4. Can one find designs among matrix groups? Section lllll 
treats this special case. It turns out that the theory is 
especially clear when one restricts attention to groups. 
The frame potential will be re-interpreted in terms of 
basic character theory. 

5. Is it possible to explicitly construct approximate unitary 
designs? In Section ITVB I we give an explicit construc- 
tion for turning mutually unbiased bases (MUBs) into 
unitary matrices which asymptotically approximate 2- 
designs. More precisely, the prescription yields a set 
of unitaries for every prime-power dimension d. These 
sets approximate 2-designs as d — > oo. Also, the car- 
dinality of such asymptotic designs is of the same order 
as the lower bound derived earlier. 

6. The set of operators B(<D d ) on <D d form a d 2 - 
dimensional Hilbert space. What is the connection be- 
tween the unitary designs in U (d) and spherical designs 
in <C d ? The relation can be made rather explicit in 
terms of the Jamiolkowski isomorphism and the frame 
potential. Both spherical and unitary designs corre- 
spond to minima of the potential - yet under different 
constraints. This statement is made precise in Section 

urn 

7. What about more general concepts such as t-designsfor 
t > 2 or substituting U(d) by other groups? We will 
discuss this general scenario in Section lVAl and present 
an example of a qubit 5-design. 
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II. GENERAL THEORY 
A. Preliminaries 

In this section we are going to derive a simple criterion for 
identifying 2-designs as well as lower bounds for the number 
of elements K they need to contain. Before stating these re- 
sults, let us shortly recall some general facts about twirling 
channels and completely positive maps which will be needed 
in the sequel. 

Let {Ug}g£G be a unitary representation of some group G 
on a Hilbert space H. The twirling channel induced by G and 
the representation U g is 

T(A) = f U g AUl dg, (5) 
Jg 

where dg stands for the Haar measure of the group G. De- 
note the projection operators onto the irreducible subspaces 
of {Ug} g by {Pi}i- For simplicity we assume that the repre- 
sentation is a direct sum of inequivalent irreducible ones (see 
Appendix I VIII Al for the general case). By Schur's Lemma 
T(A) = A if and only if A is a linear combination of the Pi's. 
Indeed, setting P( := Pi/ 'tr Pi, one easily checks that 

T{A)=Y,^{P' l A)P l . (6) 

i 

Setting H = <D d ® <C d , G = U(d) represented as U h-> 
U ® U, we arrive at the UXJ -twirling channel Tuu defined 
in Eq. (O, which has played a prominent role in quantum 
information theory. In order to identify the irreducible sub- 
spaces, define the flip operator F which acts by permuting 
the tensor factors: F|i) <8> \ j) = \j) <S> Its eigenspaces are 
the sets of symmetric and anti-symmetric vectors respectively. 
The projection operators onto these spaces will be denoted 
by P s = (1 + F)/2 and P A = (1 - F)/2. We have that 
dimPs = d{d+ l)/2 and dimi^ = d(d- l)/2. 

Moving on, we recall the well-known correspondence be- 
tween completely positive maps (cp maps) sending B(<C d ) — > 
B(€ d ) and states on € d <8> € d . Let A be such a map. Choose 

a basis {\i)}i in <C d and let |*) := Yh K) ® K) be an ( un " 
normalized) maximally entangled vector in <C d £g> <D d . The 
object 

C A := (1 ® A)|*)(*| =X)|»)0'l ® A(|*)01) (7) 

is called the C/10/ matrix of A. It is also known as the process 
matrix and the correspondence in Eq. (|7]i goes by the name of 
Jamiolkowski isomorphism. The name is justified as Am Ca 
is invertible: 

A(|i)01) = <i|iC A |j>i. (8) 

In what follows, we will write T-p for the channel induced by 
a set of unitaries T> via Eq. (O and denote the correspond- 
ing Choi matrix by Cx>- Likewise, Cjju designates the Choi 
matrix of Tjjjj. 



B. The frame potential 

The various V-quantifiers in Definition Q] make it hard to 
identify a given set of matrices as a design. Any exploration 
of this structure would thus greatly benefit from a simple cri- 
terion for the property of "being a design". Indeed, for the 
case of spherical designs such a tool is well-known (see Ref. 
Q3£| and references therein): a set of vectors {|t/>i), ■ ■ ■ > IV'-R')} 
is a spherical 2-design in <C d if and only if 

Y J \(MM\ A /K 2 = 2/{d i + d 2 ). (9) 

k,k' 

The expression on the left-hand side has been linked in Ref. 
dH] to a concept which appeared in the context of frame theory 
in an equally insightful and enjoyable paper by Benedetto and 
Fickus 13111 . The authors considered a physical model to in- 
troduce a notion of "evenly distributed" vectors: if we assume 
that K particles on the unit-sphere with respective coordinates 
\ipk) are subject to a repulsive force proportional to (V'fclV'fc') 2 ' 
then the left-hand-side of Eq. (O gives the potential of the 
configuration. Consequently, the quantity is referred to as the 
(spherical)/rame potential [46]. It turns out that 2/ {d A +d 2 ) is 
the lowest value the frame potential can possibly attain and so 
there is a one-one correspondence between global minimizers 
of the frame energy and spherical 2-designs. 

Our first result transfers this nice concept to the setting of 
unitary designs. 

Theorem 2 (Frame potential). Let T> = {Uk}k=i,...,K be a 
set of unitaries. Define the frame potential ofT> to be 

V(V)= \^ulu k >t/K 2 . (10) 

u k ,u k ,ev 

The set T> is a unitary 2-design if and only ifV (T>) — 2, which 
is a lower bound to the global minimum of the potential. 

Theorem [2] allows us to discuss the connection between 
unitary and spherical designs quite explicitly. Recall that 
the Hilbert-Schmidt inner product on B((D d ) is defined as 
(A\B)hs : = tv(A^B)/d. In the spirit of the Jamiolkowski 
map, we can establish an isomorphism between B(<C d ) as a 
vector space and <D d ® (D d . Explicitly, we map U to \vu) 
where 

v i ^=U i j /Vd. (11) 

Here, we have used the notation = (8 \y) and 
U l j = (i\U\j) for the respective matrix elements [47]. One 
checks that \vu) is a normalized maximally entangled vector 
if and only if U is unitary. Hence, we can re-phrase Theorem 
|2]as: I? is a unitary 2-design if and only if 

E \(vu k \v Uk ,)f/K 2 = 2/d\ (12) 
u k ,u k ,ev 

which is the global minimum of the spherical frame potential 
for K maximally entangled vectors. Note the close similarity 
to Eq. ©. 
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The relation between unitary designs in U(d) and spheri- 
cal designs in (D d ® <C d now becomes apparent: both corre- 
spond to minima of the frame potential, yet under different 
constraints. For spherical designs the minimum is taken in the 
set of all normalized vectors; whereas in the unitary case one 
demands that the vectors are also maximally entangled. 

Theorem |2] also facilitates numerical searches for designs 
on low dimensional spaces. Indeed, the authors have written a 
program for the MatLab computer system, which numerically 
minimizes the frame potential of a set of operators on (D 2 . If 
the set has K > 12 elements, a multitude of unitary 2-designs 
is found, while there seem to be no solutions for K < 12. 
These findings support Conjecture |4] 

Proof, (of Theorem^ Using the notation introduced in Sec- 
tion InA] let A := Ct> — Cuu- Obviously, I? is a 2-design if 
and only if ||A||| := tr|A| 2 = 0. We compute 

tr(AtA) = tviC^Cuu - C uv Cv ~ C^Guu + C^Od) 

and treat the terms in turn. To that end introduce a basis {\i}} 
in (D d <g) (D d such that the first d s := dimP s = d(d + l)/2 
vectors are symmetric and the last d a := diruP^ = d(d — 
l)/2 ones anti-symmetric with respect to F. In the formulas 
below, we will sometimes write \is) or 1^) to indicate the 
subset a given vector belongs to. Note that the vector 
introduced in Section lHAl can be written as |^>) = J^. \is) ® 
+ J2i A <S> |u>- One then finds 

Cuu = E \is)(js\®tr{\i s )(js\P s )Ps + S " A 

is ,js 

= Pa®P' a + Ps®P's- 

We have used the abbreviation S «-> A to denote the term 
which follows from the preceding one by a straight-forward 
substitution of symmetric by anti-symmetric expressions, and 
as before P' s = P s f tr P S ,P' A =P A / tr P A . Further: 

cv = X)i*>oi® (j2 u k 2 \miutr 2 /K\. 

i,j \ k ) 

Hence, 

tr(C^CW) = d s 2 tr{P s ®Ps) + d~ A 2 tv{P A ®P A ) = 2 
and 

tr(Cl v C v ) 

= K-^HPs\m) ^ (p^U® 2 \i){j\(Ut)A 

i,j \ k / 

+S ^ A 

= E tr f E U f P 's\is){is\(ut)A +S^A 

is V k ) 

= J2tr(P' s \ ls )( ls \) + S^A 

is 

- 2 = tr(C&CW). 



Lastly: 

= ^ 2 EE tr (^ 2|j)(z|( ^ t)02 u ^\( u l) m ) 

i,j k,k' 

= K-*Y,\teUtU k ,f = V{D). 

k,k' 

The claim is now immediate. □ 

For the construction of approximate unitary designs in Sec- 
tion [IVB] we record the following corollary. 

Corollary 3. Let T> be a set of unitary matrices, Cuu ond C-p 
as defined in Section Ul A\ Then 

\\C uu -Cv\\l=V(V)-2. 

C. A lower bound 

Intuitively it is clear that constructing unitary designs be- 
comes more challenging the fewer elements K one allows for 
(c.f. Theorem |6l61 l. This section is devoted to finding a lower 
bound for if as a function of the dimension d. 

What is the situation to date? Before the present paper, 
all known families of 2-designs were subgroups of the Clif- 
ford group Cd in prime-power dimensions d. In the context 
of quantum information, the Clifford group is the set of uni- 
taries mapping the set of Weyl operators (also known as: gen- 
eralized Pauli operators) to itself under conjugation l20ll . An 
introduction into this theory will be given in Section |IV CI 
where all claims made in this paragraph will be elaborated 
on. References HHEUM] use the fact that the full Clifford 
group Cd constitutes a 2-design. However, as the cardinal- 
ity of Cd grows exponentially in d (c.f. Eq. (|50l)), one might 
hope for the existence of more optimal designs. Fortunately, 
in Ref. [12] it has been realized that a particular subgroup of 
the Clifford group already possesses the 2-design property. 
The group's order scales as 0(d 5 ). What is more, the exis- 
tence of Clifford 2-designs with d 2 (d 2 — 1) = 0(d 4 ) elements 
has been established for several dimensions [12]. As will be 
explained in Section [TV CI d 2 (d 2 — 1) is in fact the smallest 
value a design based on the Clifford group can possibly have 
and that value will subsequently be referred to as the Clifford 
bound. For various reasons to be explained later, we believe 
this to be a general lower bound for the cardinality of any 
2-design, even for constructions which are not based on the 
Clifford group. 

Conjecture 4. The Clifford bound 

d 4 - d 2 (13) 
is a lower bound for the cardinality of any unitary 2-design. 

While we were not able to prove this conjecture, an estimate 
which equals the Clifford bound in leading order is established 
below. 
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Theorem 5 (Lower bound on K). A unitary 2-design in di- 
mension d has no fewer than 

d A -2d 2 + 2 (14) 

elements. 

Note that a spherical 2-design in <C d <8 <C d has at least d 4 
elements - so the slightly higher frame potential characteristic 
for unitary 2-designs might allow one to save a few elements 
as compared to spherical 2-designs. 

Proof. We follow an idea from Ref. J3]. Let \w) '■= U (8 
1\vq) for some maximally entangled vector \vq). Define a 
homogeneous polynomial p of degree 2, 2 in U by 

p(U) := (vulAl^trQvuXvulB). (15) 

Certainly, Eq. (Q3 holds and hence, as A and £> are arbitrary, 
the relation 

J2(vu\A\vu)\vu)(vu\/K (16) 
uev 

= j{vu\A\vu)\vv){vv\dU=:k(A) 

must hold and defines a channel A i— > A(A). We want to 
compute the kernel of A. 

The channel A is clearly U <8 1-covariant: 

A((f7 <8 t)A{U <8 l) f ) = {U®1) A(A) {U <8 l) f . (17) 

But because for any maximally entangled state (1 <8 V)\v) — 
(V <8 for some V', A is also 1 <E> V and hence even 

U <8 y-covariant, for all unitaries U, V. Invoking Schur's 
Lemma one concludes that A must be a mixture of projec- 
tions onto the U <8 ^-invariant subspaces of B(H). What are 
these spaces? Let us first identify the irreducible components 
of U ■ U>. The multiples of the identity (Mi for short) clearly 
form an irreducible component by themselves. Its comple- 
ment is the space of trace-less operators (M 2 ). Now, U ■ W 
must act irreducibly on M 2 , because there exists bases of mu- 
tually conjugate trace-less operators (48J. 

Surely then, Ml <£> Mi (multiples of the identity), Mi <8 
M 2 ,M 2 ®Mi (the local observables of the form 1® X, X ® 1) 
and M 2 ® M 2 are invariant under U® V ■ EH ® V'. A moment 
of thought reveals that they are irreducible (think of cyclic 
vectors). Clearly, Mi ® M\ has no non-trivial intersection 
with ker(A), while Ml (8 M 2 , M 2 ® Mi C ker(A). What 
about M2 (8 M2? Because of A's structure, either any element 
of Mi ® Mi is in the kernel or else, none is. Setting B = A* 
and evaluating Eq. (15[ we see that 

p(U) = \(vu\A\vu)\ 2 >0, (18) 

so A(A) = if and only if (u|j4|«) = for all maxi- 
mally entangled vectors |t>). To conclude that M2 (8 M2 has 
no intersection with ker A, we only need to assure the exis- 
tence of a single traceless observable X and a single maxi- 
mally entangled state \v) such that (v\X <£> X\v) ^ 0, which 
is trivially possible. Hence rank A = d A — dim ker A = 
d 4 - dim(Mi ® M 2 ) - dim(M 2 ® M x ) = rf 4 - 2(rf 2 - 1). 
But the rank of A cannot be larger than K, by Eq. ( fT6b . □ 



III. GROUP DESIGNS 

When searching for unitary designs, it might prove helpful 
to assume some additional structure in order to narrow down 
the search space and simplify the proofs. Indeed, sets of uni- 
tary matrices appear most naturally as representations of finite 
groups and (except for our numerical findings), all known de- 
signs are matrix groups. It will turn out that the concept of 
unitary designs has a very natural interpretation in terms of 
representation theory. 



A. Irreducible constituents 

We will be concerned with sets T> of unitaries which form 
a finite matrix group on <D d . It will prove convenient to con- 
ceive T> as the image of a representation U : g 1— ► U g of 
some finite group G. Recalling the notions of Section ITTAl it 
is clear that the channel Tt> is nothing but the twirling chan- 
nel associated with the representation U. By Eq. (0, Tjj will 
project onto the irreducible subspaces of this representation. 
As any operator of the form U g ® U g commutes with the flip 
operator F, we know that the symmetric and anti-symmetric 
subspaces of (D d (g> <C d will be among the invariant subspaces 
of {Ug (8 U g I g € G}. In general, these spaces are not going 
to be irreducible. We now see what makes representations U 
which induce a 2-design special: Tj> = Tjjjj (and hence T> 
is a design) if and only if the representation g 1— > U g ® U g has 
exactly two irreducible components. 

Simple as this observation may be, it must not be underes- 
timated: it allows us to understand designs from a group the- 
oretical point of view. The next section will further elaborate 
on this approach. 



B. Characters 

Let us devote one paragraph to recall some very basic no- 
tions and results from representation theory [29]. To every 
unitary representation U : g 1— > U g of a finite group, one 
associates its character ((g) = tr U g . One says that the rep- 
resentation affords £. Denote the irreducible representations 
(irreps) of G by {V^}i and their associated irreducible char- 
acters by {xi}- One introduces a scalar product between char- 
acters by setting 

(C^HGT^CteMfl)- (19) 

a 

It is a well-known and fundamental relation that the irre- 
ducible characters are ortho-normal: (xi,Xj) = The 
fact that any representation reduces to a direct sum of irreps 
means that any character can be expanded in terms of the ir- 
reducible ones and further that (C,Xi) gives the number of 
times rii the i-th irrep occurs in the decomposition of the 
representation affording Q. Finally, if Q = J^i n iXi^ then 

IICII 2 = (CO 
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Now, let T>, G, U be as in Section MI A I We compute the 
frame potential of T>: 

V(V) = ]T|tr[V v | 4 /|G| 2 (20) 

g,g' 

= £|trf///|G| 
g 

= ( tr U 9 ® u g)(& U s ® U g)/\ G \ 
g 

= (Cc/<2) ^ Ct/C2) ) = ||C[/(2)|| 2 , 

where Cu<. 2 ) (g) = tr(U g ® U g ) — Cu(g) 2 is the character of 
the representation : g > U g <£> U g . In other words, the 
frame potential of a group design is the squared norm of the 
character of U^ 2 \ 

We can now rederive Theorem [2] By this section's first 
paragraph, ||C[/( 2 )II 2 = 2~2i n h which equals 2 if and only 
if U^ 1 has exactly two irreducible components. This in turn 
is equivalent to U inducing a 2-design, as has been shown in 
Section UlIAl Note how much the group structure simplified 
the proof. 

It seems remarkable that the frame potential offers a very 
natural interpretation in terms of two completely unrelated 
structures: from the point of view of frame theory, it is a 
purely geometrically motivated measure for the "eveness" of a 
distribution. In terms of group representation theory, it seem- 
lessly takes on an algebraic role. 



C. General results and properties 

Consider a group design D — {U g \g £ G}. The cen- 
ter Z(U) of U are the elements of T> which commute with 
any U^. By Schur's Lemma, if U is irreducible, we have that 
U g € 2.(1]) Ug oc 1. Now choose one representative of 
each coset T>/Z(U) and assemble these unitaries in a set T>' 
(D' is called a transversal of T>/Z(U)). Using Eq. (0, one 
sees that D' is a 2-design of cardinality |2?|/|Z(£/)|. From 
now on, we will restrict attention to such reduced sets. Con- 
sequently, for any representation U of G, we will define T>jj 
to be a transversal of {U g | g S G} / Z(U) and refer to T)jj as 
the group design induced by U . 

Using this definition, let us collect and extend the results on 
group designs in the following theorem. 

Theorem 6 (Group designs). Let G be a finite group and U a 
unitary representation of G on <C d affording the character 
The following are equivalent: 

1. The set T>jj is a 2-design. 

2. The representation U^ : 5 n U g <g> U g has no more 
than two irreducible components. 

3. It holds that \\Cum 1 1 2 = 2 - 



4. The characters 

xs(g) ■■= «5) 2 + x(s 2 ))A 
XA(g) := (x(g) 2 ~x(g 2 ))/2 



are irreducible. 
Further: 

5. The cardinality K — \T)jj\ is a multiple of d and 
l/2d(d±l). 

6. Let H be a finite group represented on <C d by V. If 
{V h I h e H} D {U g I g G G} then T{V V ) < T{V V ). 

7. The frame potential of a matrix group is an integer. 

8. A necessary condition for T>u to be a 2-design is that U 
is irreducible. 

9. For d > 2, there are no real-valued group representa- 
tions which form a 2-design. 

Statement [5] can be used in conjunction with Section HTCl 
to derive bounds on K. For example, for d = 2, it holds 
that K > 10 by Theorem[5] But we now know that K must 
be divisible by 2 and 3, so that K > 12. As unitary group 
designs of order 12 in dimension 2 do indeed exist, we know 
that the bound is tight. Unfortunately, this is the only case 
where we can make lower and upper bounds match. Note 
also, that 12 is the value predicted by the Clifford bound for 
d = 2, supporting our conjecture. 

The|6}th point says that "supergroups have lower frame po- 
tential than their subgroups". Again, it is clear that construct- 
ing designs is easier, the more elements one allows for. In gen- 
eral, however, just adding further unitaries to an "almost de- 
sign" is not going to improve the potential. For group-designs 
the situation is different, as we now know. 

Lastly, statement [7] says that the frame potential of matrix 
group is "quantized". In that sense, there are no "approximate 
group designs". 

Proof. The equivalence 1. 2. 3. has been established 
in the discussion preceding the theorem. Claim 4. is equiv- 
alent to 2., as x 2 — Xs + Xa- The fifth statement follows 
from a well-known theorem in representation theory (see Ref. 
l29ll ). Point 6. holds true as the number of irreducible com- 
ponents cannot decrease when passing from a subgroup to a 
supergroup. Claims 7. and 8. should be obvious. Lastly, 9. is 
valid because for real (u 

1 = (CuXu) = (CuCu, 1g) = (C[/< 2 > I 1g?)j 

where 1^ : g *—>■ 1 is the trivial representation. Hence, 1q 
is a one-dimensional irreducible component of U^ 2 \ But for 
d > 2 we have that d s ,d a ^ 1. □ 
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D. Harvesting character tables 

The results of Section IIIIBI enable us to identify designs 
by just looking at character tables of finite groups. Such ta- 
bles have been the subject to intensive research and are digi- 
tally available. We have employed the freely available GAP 
computer system 113211 to search the GAP Character Table Li- 
brary version 1.1 IB3I1 for unitary designs. Some findings are 
compiled in Table U For each dimension d in which a uni- 
tary design has been found, one example is included in the 
table. To access the listed character tables, pass the name to 
CharacterTableFromLibrary ( ) . The column 'Trred. 
character no." gives the position of the design within the list of 
irreducible characters returned by the Irr ( ) -function. E.g., 
the dialog 

gap> t : =CharacterTableFromLibrary ( " J 4 " ) ; ; 
gap> chr :=Irr (t) [2];; 
gap> Degree (chr); 
1333 

gap> Norm (chr*chr) ; 
2 

gap> 

confirms that the last item in Table U does indeed belong to a 
unitary group design in dimension 1333. 

IV. PHASE SPACE TECHNIQUES 

The title phase space techniques refers to any method em- 
ploying the realted concepts of Weyl operators (also known as 
generalized Pauli operators), stabilizer states and the Clifford 
group. These structures have played a central role in the the- 
ory of both spherical and group designs J2, H, Gil EH H • As 
the following paragraphs require some rather technical prepa- 
rations, we state a summary of the results at this point. 

In Section IIV Bl asymptotic unitary designs will be con- 
structed. By this, we understand a family of sets of uni- 
taries T>d, such that the matrices in T>d are d-dimensional and 
lim^oo ViJ^d) = 2. The intuition behind the construction 
is as follows: in Section iHBl we discussed the relation be- 
tween the frame potential of operators on (D rf and vectors in 
<C d . Hence it is natural to ask whether one can exploit this 
relation to turn spherical designs into unitary ones. Obviously, 
in order to obtain unitary matrices, we must require the vec- 
tors in the spherical design to be maximally entangled. Recall 
that a maximal set of MUBs is a 2-design and, moreover, that 
such sets can be chosen to consist of stabilizer states ifioll . For 
bi-partite systems, where each party has prime dimension, it 
is known that stabilizer states are either maximally entangled 
or not entangled at all l30ll . It is thus reasonable to assume 
that among the elements of a maximal set of MUBs, there are 
"enough" maximally entangled ones to yield a set of unitaries 
with a low frame potential. Fortunately, this intuition turns 
out to be true and we will find sets of 0(d 4 ) unitaries in di- 
mension d = p n , which approximate a 2-design as d — > oo. 

Secondly, in Section |TV CI we will revisit the technique of 
Clifford twirling. Our main contribution to the theory will be 



a systematical reassessment of what is already known. Indeed, 
reading the literature, one gets the impression that some con- 
fusion has arisen due to the fact that several distinct Clifford 
groups exist. The one used in Refs. HI 0. EH HH] is differ- 
ent from the one in Ref. 11211 . Going on, we will review a 
construction by Chau, which meets the Clifford bound in di- 
mensions 2, 3, 5, 7, 11, and outline a way for circumventing a 
no-go theorem which asserts that for any other dimension the 
bound cannot be met. In particular, for d — 9, we present a 
subgroup of the Clifford group which is a 2-design of smaller 
cardinality than Ref. [ 12] seems to suggest is possible. 

Before presenting these results in detail, the reader must en- 
dure the tour-de-force of technical preparations given in Sec- 
tion llV Al It is a peculiarity of the theory to be presented that 
it works much more smoothly in odd dimensions d than in 
even ones. While it can be checked that all results in the next 
section also hold for the qubit case, the proofs are given only 
for the case of odd d. 



A. Introduction 

This section contains a very brief outline of the general the- 
ory. See Ref. l2~5l l26ll and references therein for a more de- 
tailed exposition. 

1. Weyl operators, the Jacobi group & the Clifford group 

Let us first gather some well-known facts on finite fields 
l28tl . If p is prime and r a positive integer, F pm denotes the 
unique finite field of order p m . The simplest case occurs for 
m = 1, when F p ~ 7L V , i.e., the set of integers modulo p. 
Now set d — p m and choose an r € N. Out of the base 
field B :— Fd, one can obtain the fields F^ by means of 
a field extension. Extension fields contain the base field as 
a subset. The extension field possesses the structure of an r- 
dimensional vector space over the base field. A set of elements 
of F is a basis if it spans the entire field under addition and 
^-multiplication. The operation 

Tr F/s / = ][><** 

fc=0 

takes on values in the base field and is _B-linear. Therefore, 

(f,9) !-> Tr F/B (fg) 

defines a B-bilinear form. For any basis {hi], there exists 
a dual basis {b 1 } fulfilling the relation TrpmQfbj) = Sij 
(we do not use Einstein's summation convention). Clearly, if 
/ € F can be expanded as / = J^. with coefficients 
/' S B, then duality implies that /' = Tr(/6 4 ). 

We will work in the d := p m -dimensional Hilbert space 
H ~ <C d spanned by the vectors {\a) \ a S F^}. Define a 
character of F d by Xd(a) ■= exp(i^ Tr Fpm / Fp (a)). The 
relations 

Xd(q)\x) = \x + q), z d (p)\x) =Xd(px)\x) (21) 
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define the shift and boost operators respectively. The Weyl 
operators (also known as generalized Pauli operators) in di- 
mension d are given by 



wd{p, q) = Xd(-2 1 pq) z d (p)x d (q), 



(22) 



for p, q € ~Fd- The phase factors in Eq. d22t have been in- 
cluded to clean up some later formulas. The phase space V 
is defined as V := x F^. We introduce the standard sym- 
plectic inner product on V by 



P 



pq -qp 



(23) 



For elements a = (p,q) T of V, we set Wd(a) ■— Wd(p,q). 
Denote by Wd '■= {wd(a) | a E V} the collection of all Weyl 
operators. The commutation relations 



Wd(a)wd(b) = Xd([a,b])wd(b)wd{a), 



(24) 



can be checked to hold. 

Let S be a symplectic 2x2 matrix with entries in Wd- There 
exists a unitary operator fid(S) defined via 



P>d(S) w d (a) Hd(S) 1 = w d (Sa) 



(25) 



for all a E V. We will call fid(S) the metaplectic representa- 
tion of S. Up to phase factors, the set of unitaries of the form 
Hd(S) uid(a) constitute a group, which will be referred to as 
the Jacobi group 

The preceding definition have been made with a single 
d-dimensional particle in mind. We now consider the sit- 
uation of n particles, each having d — p m levels. The 
HUbert space becomes <C d spanned by {\a) | a E F^}. Let 
p = (pi, . . . ,p n ),q = (qi, ■ ■ ■ , q n )- We define the Weyl oper- 
ators as 



Wd,n(p, q) = Wd(pi,qi) • • • Wd(p n , q n )- 



(26) 



In this case, the phase space is set to be V = FJJ x F^ and 
Eqs. ( 1231241 1 continue to make sense if we perceive products 
between elements of p, q E as a canonical scalar product: 
pq = Y^? PiQi- I n complete analogy to the n = 1 case, one 
finds that for any symplectic 2n x 2n matrix S with entries in 
Wd, there exists an operator fid,n(S) such that 



Pd,n(S) W d ,n{a) Hd,n(Sy = W d ,n(Sa) 



(27) 



holds for all a E V. Denote the set of Weyl operators ac- 
cording to Eq. d26"l ) by Wd, n an d the Jacobi group spanned by 
{^(iKii^Ks by Jd, n - 

For a HUbert space of prime-power dimension p s , we can 
now construct an entire family of different Weyl operators and 
Jacobi groups. Indeed, for any n, m such that nm = s, the 
Weyl operators w pm n are p s dimensional. Prominent choices 
include n = lj_m = s (used in Ref. |[l2tl ) and n = s, m = 1 
(used in Refs. ylEHIZlI]). It will turn out that all definitions of 
the Weyl operators coincide, while the various Jacobi groups 
differ. Proposition|7]makes these remarks precise. In order to 
state it, we need one final definition: the Clifford group C P) „ 
is the set of unitaries mapping the set W pn onto itself under 
conjugation. This definition reflects the general use of word 
Clifford group in quantum information theory l23ll . 



Proposition 7. Let p s be a power of a prime. Let n < n' and 
to, to' be such that ran — m'n' = s. Then 



W P ™,n = >Vv 
The inclusion in Eq. ( 1291 ) is proper and J v . s = C p ., 



(28) 
(29) 



For the construction in Section IIV Bl it will be necessary 
to understand Eq. d28| i in more detail. Indeed, it has been 
realized before I25l l26l 12711 that the Weyl operators in W p ^ i 
can be written as tensor products of those in W p , n . In what 
follows, we will refine this picture. 

Let d = p m be a power of a prime, let B = F^. Let 
F = Frfn be an extension field of B. In F, choose a basis 
{bi}i—i...n over B. Denote the dual basis by {&*};. Having 
general relativity conventions in mind, we will adopt the fol- 
lowing notation: for an element / E F, we denote its expan- 
sion coefficients with respect to bi by /' and the coefficients 
for the dual bases by /*: 



(30) 



The Weyl operators in W d ,i act on H = C d " ~ (<C d f 
where we choose the isomorphism to be implemented by 



\q) = tfh + ■■■ + q n b n ) ^ k 1 ) 



\q n ). (3D 



Lemma 8 (Factoring Weyl operators). Using the notions in- 
troduced above, the Weyl operators in Wd™ ,i factor as 



Wd"(p, q) = WdiPi^ 1 ) 



w d (Pn,q n )- (32) 



Proof. Denote the common prime field F p of B and F as P. 
It is well-known 12811 that Tr B / P o Tr F / B = Ti F / P . Hence 

XFipq) = XBi^Pjq'^F/Bihb 3 )) = Y[xB(Piq 1 )- 



Similarly, 

i j i 

= (g)^(g l )|x l ), 

i 

i j i 3 

= ^)z B { Pt )\x l ). 



Using Eq. d22l i. the claim follows. 



□ 



Proof, (of Proposition^ Eq. d28l follows from the previous 
lemma. Saying that J p r = C Ptr is just rephrasing the defini- 
tion of the Clifford group. For Eq. d29l ) the reader is deferred 
to Refs. HI US. □ 
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2. Stabilizer states 

Using the commutation relations Eq. $2M it is immediate 
that to Weyl operators w(a),w(b) commute if and only if 
[a, b] = 0. Now consider the image of an entire subspace 
M of V under w: 

w(M) = {w(m)\m £ M}. 

The latter set consists of commuting operators if for all 
m,i,TO2 € M, the symplectic inner product vanishes: 
\rn\,m,2\ = 0. If that condition is fulfilled, w(M) is called 
a stabilizer group. Consider the operator 

Pm := w(m)/\M\. (33) 

One checks that 

trp M = ^trto(m)/|M| (34) 
m 

= ^ dS ^o/\M\=d/\M\. 

tn 

and, using the fact that M is a linear space, 

PMPM = \M\~ 2 ^ W ( m + m> ) ( 35 ) 

m , m ' 

= |Af| _1 5^«>(m) =pm- 

m 

Hence, if \M\ = d, then pm = \4 1 m){'4 ) m\ is a rank-one 
projector and \ipni) is called the stabilizer state associated 
with M. The preceding definition can be extended: choose 
a character £ of M (i.e., a function M — ► (D such that 

£(mi + m.2) = C( m i)C( TO 2)) an d set 

Pw,C : =2C(m)«;(m)/|M|. (36) 

m 

The calculations Eqns. ( 134135b can be repeated and one finds 
that also pM,C projects onto a vector, which will be denoted 

by \iPm,c)- 

3. Mutually unbiased bases 

We will recall a well-known construction for MUBs floll . 
Once again, let F = F p m be a finite field and V — F 2 the 
associated phase space. Let 

Va=(f), M a = {\v a \aeF}. (37) 

Clearly, M a is a one-dimensional subspace of V and hence of 
cardinality \M a \ = \F\ = d. Because the symplectic form is 
anti-symmetric [a, b] — —[b, a] it holds for Xv a , \'v a £ M a 
that [Xv a , X'v a ] — \A'[v a , v a ] — and hence the spaces M a 
fulfill the requirements of the last section and define stabilizer 



states. Further, setQ a '(Xv a ) := &(A) : = exp(i^ Tr(6A)). 

Each is easily seen to be a character of M a . Now define 
the stabilizer states 

B£ o) £ Ci(A)«>(Aw«). (38) 

A 

We claim that these states constitute a set of MUBs. Indeed 

trB^B^ - d- 2 ^C b (A)a'(A')tr W (A« a + AX) 
XX' 

= ^ 1 X)C6(A)Cv(A / )«a«.,ax ^ 

XX' 

Now if a = a' then Eq. d39b reduces to 

tr Bi a) B^ = d-^exp^— TV(A(6-6'))) 

A P 

= <S 6 ,f/ (40) 

while for a ^ a' we use the property of the finite plane F 2 
that the two lines M a and M a i intersect exactly at A = to 
conclude 

trB^B^ = d- 1 C 6 (0)Ch'(0) = d- 1 . (41) 

Hence, for a fixed a, the set {B^}b forms a basis and all 
d bases corresponding to different values of a are mutually 
unbiased. The computational basis corresponds to the set 
Moo = {A(1,0) T |A S F}. Repeating the reasoning em- 
ployed above, one finds that it is unbiased with respect to all 
the other ones; hence, we have constructed a maximal set of 
d + 1 MUBs. 

B. Asymptotic designs from MUBs 

This section revolves around the following definition. 

Definition 9 (Asymptotic 2-designs). Let I C N be an index 
set. A family of sets of unitaries T>^, d £ I is an asymptotic 
2-design if the matrices in T>d are d-dimensional and 

lim V(V d ) = 2. (42) 

d — >OG 

A priori, it is not clear that this definition has any physi- 
cal relevance. After all, it is conceivable that, even though 
the frame potential of T) d converges, the XVtwirling chan- 
nels Tj) d do not become close to the {/{/-twirling channel in 
any sensible metric. Indeed, the question of whether asymp- 
totic designs are "almost as good" as strict ones cannot be 
answered in general, but depends on the application one has 
in mind. One particular aspect of this question will be illumi- 
nated in Lemma [TO] We will show that the series of twirling 
channels T-p d does converge to Tuu in Z? plo -norm. The lat- 
ter norm has been defined in Ref. |35], a well-readable ac- 
count of the merits and perils of different metrics for quan- 
tum channels. Specifically, let A and A' be channels with 
respective Choi matrices C,C. If A = C — C , then 
D pm (A, A') := d _1 tr|A|. Some physical interpretations of 
Dp ro -convergence are listed in Ref. [35]. 
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Lemma 10. Let T>d be an asymptotic 2-design. Then the T> e i- 
twirling channels Tx> d converge to Tyu in D pm -norm. 

Proof. Let C-p d , Cuu be the usual Choi matrices, let A = 
Cx> d — Cjju- As quantum channels preserve Hermiticity, the 
Choi matrices and hence A are Hermitian. Let {Si} be the set 
of eigenvalues of A. By CorollaryO we know that tr |A| 2 — > 
0. Hence, for all i, |<J,-| < 1 eventually and therefore, for large 
enough d: 

D pro (A d ,A a[/ ) = d- l Y,\^\ <d- 1 Y J \^\ 2 

i i 

= d-^IAlH^O. 

□ 

The technically non-trivial part of this section is contained 
in the next theorem. 

Theorem 11 (Mutually unbiased bases). Let d — p m be a 
power of a prime. Then in H — € d <£> € d exist d 2 + 1 MUBs 
of which d 2 — d bases are maximally entangled and d+1 bases 
factor. 

Proof, (of Theore mJIB Le t B — ¥4, F — F d 2. Using the 
notation of Section ITV A II we assume that a basis for F over 
B has been chosen. For a,b £ F, let B^ be a projection onto 
a stabilizer state in fi, as defined in Section llV A 31 Hence 

B W = d -2 J- (l a) (p,q)w F (p >q ) 

(p,q)£M a 

= dr 2 ^2 Ci a \p^) w B(pi,q 1 )®WB(p2,q 2 )- 

{p,q)eM a 

Defining N a = {(p, q) e M a \p 2 =q 2 = 0}, we get 

tr 2 ^ a) = cT 1 Y, ^\ Pl q)w B (p 1 ,q 1 ). (43) 

(P,«)GA f a 

Clearly, N a is a £>-vector space, so it has cardinality |iV a | = 
d n for some n. If n = 0, then tr 2 B^f 1 equals 1 1; so the state 
was maximally entangled. If n = 1, Eq. (l43l is of the form 
of Eq. ( 1361 ) which makes tr 2 B^ a pure state on <C d . Further, 
n = 2 would imply |JV | = |Af | and hence B^ = pi ® 1 2 

for some density operator pi, which is impossible as Bi is 
pure. Lastly, n > 2 \N a \ > \M a \, which is absurd. Hence 
any vector in the standard set of MUBs is either a product or 
else maximally entangled. 

Now, (aX, A) € N a if and only if Trp/ B (aAe 2 ) = 
Ti"F/s(Ae 2 ) = A = Xie 1 A aX = (aA) 1 ei. Assume 
that 

a = b% (44) 

e 1 

for some b € B. Then A = Aie 1 => Xa = Xibei and hence 
\N a \ = d. Conversely, assume that |iV a | = d. Then for all 
A = Aie 1 we must have that Xa = be\ for some b G B. 
Solving for a shows that Eq. d44t must hold. Hence among 
the d 2 bases associated with the sets M a , there are exactly 
\B\ = d factoring ones. Taking the computational basis into 
account, the assertion becomes immediate. □ 



The validity of Theorem[TT]implies the existence of asymp- 
totic 2-designs. 

Corollary 12 (Existence of asymptotic 2-designs). Let I be 
the set of prime-power integers. Then there exists an asymp- 
totic 2-design T>d, for d £ /. 

Proof. We compute the frame potential of the d 2 (d 2 — d) = 
d 4 — d 3 unitaries T>d which can be constructed via Eq. ( fTTT ) 
from the maximally entangled MUB vectors of TheoremfTTl 

v{v d ) = d^EK^I^)! 4 /^ 2 

a, a' b.b' 

= d 4 (l + {d 2 -d-l)d- 2 ) /K 

2d 4 -(d-^+d- 2 ) 
= dT^Jz > 2 

□ 



C. Clifford designs 

Let us review the technique employed in Ref. |Q1 0, fl3ll to 
construct a 2-design. The construction proceeds in two steps. 
First one realizes that twirling an operator p by Weyl matri- 
ces reduces p to its "Weyl-diagonal" components (see below). 
Secondly, twirling the resulting operator using the metaplectic 
unitaries p(S) "evens out" the coefficients to yield a U <£> U- 
invariant state. 

Denote by Tyy the Weyl twirl channel ffl [l3ll : 

T w (p) = d~ 2 Y w ( a ) ®w{a) pw(a) t <g) w(a) f . (45) 

Expanding p = >~2 b Pb.v w{b) ® w(b'), we compute: 
Tw(p) 

= dr 2 pbfi w(a)w(b)w(a)^ <8 w(a)w(b')w(a)^ 

a,b,b' 

= d- 2 Y Pb,v w(b) <8 w{b') Y X(~[b + b', a}) 

6,6' a 

= YPb,-bw(b)®w(-b). (46) 

6 

The final transformation follows from the fact that x([6, • ]) is 
a non-trivial character of V for any 6^0. 
The Clifford twirl channel Ill rOl is given by 

T c (p) (47) 

= I^f ? I~ 1 £ [u®u) P {u®u)i 

p 

= I^f-.|- 1 E (»(S)®p(S))T w (p)(p(S)®p(S)?, 

SeSp(p,n) 

where Sp(p, n) denotes the group of symplectic matrices on 
V = Fp". Using the notation and results of Eq. d46l l. one 
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concludes: 



V. MISCELLANEOUS TOPICS 



Tc(p) 

= \Sp{p,n)\- 1 Y,Pb,-b ™(Sb)®w(~Sb) 
bev ses P (p,n) 

= al + (3 ^ w(b)®w(-b), 



A. Higher orders and general groups 



o^bev 



where the constants are given by 



po,o, = Y Pb - b - 
o^bev 



(48) 



(49) 



Eq. (08]) follows because the symplectic group acts transitively 
on := {a 6 V | a =^ 0} and hence maps every element of 
V" equally often to every other element. It is evident that Tq 
projects onto exactly two subspaces and is hence equal to Tjju 
by SectionHHI 

Now let G be some subgroup of Sp(p, n) such that G 
acts transitively on VK From the above argument it is clear 
that T^g) ° Tw = Tjju and hence that {w(v)(i(S) \ v € 
V, S 6 G} is a unitary 2-design. An obvious choice is to set 
G = Sp(p™, 1) which is the basis of Ref. lHal . The advantage 
of going from the multi-particle picture to the single-particle 
picture is an exponential reduction of the cardinality of the 
design: 



n-l 



Sp(p,n)| = p n Y[(p 2{n -' l) - 1) 



(50) 



i=0 



Sp(p", 1) 



0{p 2n2+n ) = O{d 2{log p d)2+ ' los p d ), 



p n (p 2n - 1) = 0(p 6n ) = 0{d A ) (51) 



(see Ref. Bol for a derivation). What possibilities are there 
to further improve the cardinality? Clearly, any group act- 
ing transitively on must have order k |V"| = k(d 2 — 1) 
for some integer k. The smallest value is k = 1 and hence 
d 2 (d 2 — 1) gives a lower bound to the number of elements 
a Clifford design arising from such a construction can have 
(c.f. Conjecture |4j. In Ref. lfl2tl Chau showed that for 
d = 2, 3, 5, 7, 11, such minimal subgroups do exist. He goes 
on to rule out the existence of any subgroup G of Sp(p™, 1) 
which has cardinality a multiple of V*. Hence no reduction 
below d 2 \ Sp(d, 1)| = d 5 — d 3 seems to be possible in gen- 
eral. 

However, the argument leaves open the possibility of find- 
ing subgroups G of Sp(p, n) which act transitively on the 
non-zero elements of the vector space and are smaller than 
| Sp(p n , 1)|. While we do not know if such groups exist in 
general, we know of one example for d = 9. In Table |ll] 
we list the generators of a transitively acting subgroup G of 
Sp(3, 2) of order 160 = 2|F tt |. It yields a Clifford 2-design of 
cardinality 2(d 4 — d 2 ) — 12, 960, where the design induced 
by Sp(9, 1) has 58, 230 elements and the one associated with 
Sp(3, 2) consists of 4,199,040 unitaries. All claims made 
about the generators can easily be tested by a computer alge- 
bra system. 



We saw in Section UTAl that the effect of the twirling chan- 
nel T induced by some group G and a corresponding unitary 
representation U g depends only on the the decomposition of 
Ug into irreducible constituents (see Appendix IVIH A| for a 
version of Eq. (|6]l valid for general representations). 

Based on this observation, the notion of a group design can 
easily be generalized: 

Definition 13 (General group designs). Let U be some group 
of unitary matrices. Then a finite subgroup T>ofU is a U- 
group design of order t if the t-th tensor power of T> has 
the same number of irreducible constituents as the t-th ten- 
sor power of hi. 

The most natural setting for applying the above definition 
is given by unitary designs U = U(d) of higher order t > 2. 
How many irreducible constituents do we expect the represen- 
tation U i—* f7®* to decompose into? The following lemma 
answers this question by giving the frame potential for uni- 
tary t-designs, at least in two special cases. 

Lemma 14. The frame potential of a unitary t-design in di- 
mension d is given by 



[n/2j /n\(n-2i+l) s 



El=o ( 



i\(n-i+l) 



for d>t, 



for d = 2. 



Once again, we can search the GAP library for examples. 
Table [ill] gives an example of a matrix group G in d = 2, 
whose 5th tensor power decomposes into 42 irreps, the re- 
quired value for a 5-design. As a matter of fact, G is very close 
to being a 6-design: its 6th tensor power has 133 irreducible 
components, whereas the full unitary group U (2) decomposes 
into only 132 irreps. Note that the matrices in TablelTnlare not 
unitary in the standard basis. However, as is well-known i24Tl . 
any representation of a finite group is equivalent to a unitary 
one: one can easily construct a similarity transformation map- 
ping the given matrices to unitaries. 

Proof. ( of Lemma I7?l) The following facts are well-known 
rf24Tl : (i) there is a one-one correspondence between irre- 
ducible components of the t-th tensor power of U(d) and 
young frames T partitioning the integer t into no more than d 
parts; (ii) the multiplicity of the irrep belonging to a specific 
frame T is given by the dimension djr of the corresponding 
irrep of St- 

If d > t, the restriction "no more than d parts" becomes 
irrelevant. Using the results of Section ITlIBI we find that the 
frame potential of a i-design is 



E 



dj, = \S t \=tl, 



where the sum is over all young frames with t boxes. The 
second claim follows in a similar fashion, using well-known 
formulas for dp (which can be found, e.g., in Ref. [24]). □ 
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B. Energy-preserving operations as in linear optics 

We will briefly comment on a specific representation of the 
unitary group, which plays a prominent role in linear optics. 
General references for the introductory paragraphs are Refs. 
lfl6l[37ll38i l39ll . The most central mathematical object in the 
description of physical systems of d bosonic modes are the 
creation and anihilation operators a,k, at. Recall that the set of 
unitaries on L 2 (H) which keep the vector space spanned by 
the etfc , a\ invariant under conjugation form a projective repre- 
sentation of the real symplectic group Sp(2d). The represen- 
tation is referred to as the metaplectic representation. Bosonic 
systems are often thought about in terms of their phase space 
description, where these metaplectic operations appear in an 
especially natural way. 

The maximal compact subgroup of Sp(2d) is given by the 
intersection of Sp(2d) with the set of orthogonal transforma- 
tions SpO(2d) := Sp(2d) n 0(2d). Physically, the elements 
of SpO(2d) correspond to energy-preserving or passive op- 
erations. These operations are exactly the ones which are 
easily accessible in the laboratory using passive linear optical 
elements (phase shifts and beam splitters, but no squeezers, 
which correspond to elements in Sp(2d) which are not con- 
tained in 0(2d)). It is a well-known fact that SpO(2d) ~ 
U(d), which might trigger some hope that our theory could be 
applicable to these systems. However, the metaplectic repre- 
sentation is infinite-dimensional and in this work we did not 
develop the means to cope with such representations. 

In an indirect approach, we can, however, nevertheless ex- 
ploit the developed formalism: Fortunately, important prop- 
erties of bosonic quantum states can be described entirely 
in terms of objects on the 2d-dimensional phase space. In- 
deed, define the canonical coordinates or quadrature opera- 
tors r = (xi, . . .,xa,pi, . . . ,pd) by 

x k := [a k + 4)2-V 2 , p k := i{a\ - a^ 1 ' 2 . (52) 

A much-studied object in particular in quantum optics are the 
various second moments of the quadrature operators with re- 
spect to a given state. Assuming that all first moments vanish, 
the second moments can be conveniently assembled in a real 
symmetric 2d x 2d covariance matrix 7, defined as 

1 \ = 2 (tr(r fc ri) + tr(r,r fc )) • (53) 

An interaction process preserving the energy would then give 
rise to a map 

7 1 — ^ S-yS T -. 7', (54) 

where S € SpO(2d). The following proposition assures that 
unitary designs can be used to tackle problems in this context. 

Proposition 15 (Averaging over passive operations). Let V C 
U(d) be a unitary group design of order t. The image ofT) in 
SpO(2d) under the usual isomorphism is then a SpO(2d)- 
group design of order t. 

A setting that can be studied using designs in this way is the 
following: consider a system of d interacting bosons. Hav- 
ing a "microcanonical ensemble" in mind, we might be inter- 
ested in the expected value of various quantities after random 



energy-preserving interactions have been applied. Haar aver- 
ages over SpO(d) would constitute a sensible model for such 
a system (see also Ref. lfl5ll ). 

Let us give a concrete example. Take { \i) }, as a basis of the 
vector complex space in which the complex moment matrices 
are defined. It is straightforward to see that the mean energy 
of the first mode is given by E = (l|Sl7$lt|l). The quantity 

AE = [ (l\(U®U)n^(U®Uy\l) 2 dU (55) 

JU(n) 

- ( [ (i\(u®u)n^(u®u)^\i)du) . 

\JU(n) ) 

gives the expected energy fluctuations of the first mode and is 
directly amenable to evaluation using unitary 2-designs. 

Proof ( of Proposition [771) . The usual isomorphism between 
elements U E U(d) and elements S 6 SpO(2d) can be stated 
as follows: If U = X + iY, then S is given by SlU 

S(U) =( K X Y x)- (56) 

Setting 

n -(l"i)/<>* < 57 » 

one finds easily that 

nSiU)^- 1 = U © U. (58) 

As T> is certainly a group i-design if T> is, the statement fol- 
lows. □ 



C. Random entanglement 

Originally posed by Lubkin and later popularized by Page, 
the following questions has a long history: what is the av- 
erage entanglement of a composite system in a pure state 
B15[ kll l42l 14311 . One motivation for studying such problems 
is to justify the ad hoc "rule of minimal prejudice" employed 
in statistical physics, which states that among the ensembles 
compatible with macroscopic observables the one maximizing 
the entropy is realized in nature. The problem was originally 
stated in terms of the entropy of entanglement of subsystem 
A: S(pa) — — t?(pA log Pa)> but various other measures, for 
example the purity tr(p^) = °f PA can be used. The 

latter quantity has the advantage that its expectation value 

/ || tr s [Uptf] || \d\J (59) 

JU(d) 

is a Haar integral of a second-order polynomial and can thus 
be directly evaluated by averaging over a 2-design (above, p 
projects onto an arbitrary 2-system state vector \tjj})- 

Based on this observation, we obtain a simple answer to a 
special case of this problem as a corollary to TheoremfTTKc.f. 
Ref. 114411 for a much more general, but much longer proof). 
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Corollary 16 (Average entanglement). Let d be the power of 
a prime. The average entanglement of pure states on <D d <g) <C d , 
as measured by the purity, is 2d/ (d 2 + 1). 

Proof. We choose p = |0)(0| <8> |0)(0| and average over the 
Clifford group Jd 2 i- The image of |0) <g> |0) under the action 
of J(p,i constitutes of the bases of Theorem [TTI The purity 
of a reduced density matrix of a product state equals 1, for a 
maximally entangled state it is d^ 1 . We can thus compute the 
average by counting: 

(d 2 - d)d- 1 + (d + 1) _ 2d 
d 2 + 1 ~ d 2 + 1 

as claimed. □ 

Note that, while the question of finding the expected entan- 
glement of a pure state is stated in terms of analysis, we could 
answer this special case by purely combinatorial means. 

VI. SUMMARY 

In this work, we presented a first systematic analysis of the 
mathematical structure of unitary designs. We pointed out 
a connection to group representation theory, gave bounds on 
the number of elements of a design, made the relationship to 
spherical designs explicit, and used this connection to con- 
struct approximate unitary designs. Foremost, the pivotal con- 
cept of a frame potential has been explored. Intriguingly, the 
latter quantity appears very naturally in different, seemingly 
unrelated areas. 

While much insight into the structure of unitary designs has 
been gained, many questions remain unresolved as interest- 
ing open problems: finding a systematic way for construct- 
ing designs for any choice of parameters t, d or improving the 
bounds for their cardinalities, to name just two. 
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VIII. APPENDIX 

A. General twirling channels 

Eq. © gives an explicit formula for a channel twirling over 
an irreducible representation U g . In this section, we state the 
relation for the general case. So assume that g i— > U g de- 
composes into a set of irreps f/W, which have dimension di 
and occur with multiplicity n.j respectively. The underlying 
Hilbert space H then decomposes as 

W = 0Hi®C ni , (60) 

where Hi — <C di (see, e.g., Ref. IT24I] ). The representation U g 
itself can be written as 

U 9 =1>2 U 9 1) ® W (61) 

i 

Now set Pi = \<n i l n .. The twirling channel becomes 

T A (p) = J2^n i (Pip)Pi, (62) 

i 

as can be checked without difficulty. 
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TABLE I: Some group designs found by the GAP system. The group name and character number refer to the names used by the "ctllib" 
package (H]. 



2 
3 
4 
5 
6 
8 
9 
10 
11 
12 
13 
14 
18 
21 
26 
28 
41 
45 
342 
1333 



K 



12 

72 
1920 
25920 
40320 
20160 
19440 
190080 
13685760 
448345497600 
4585351680 
87360 
50232960 
9196830720 
17971200 
145926144000 
65784756654489600 
10200960 
460815505920 
86775571046077562880 



K/(d 4 -d 2 ) 



1 
1 

8 

43.2 
32 
5 
3 

19.2 
942. 
21772800 
161501. 
2.29 
480 
47397. 
39. 
237714. 
23294225607. 

2.49 
34. 

27483822. 



Group 



3.3 



'2: (3x2A4) 
2"3.L3 (2) 
2 .HSM10 
2"6:U4 (2) 

. L3 (4) .2_1 
4_1 .L3 (4) 

1 (1 + 4) : 2S5 

2. M12.2 
6xU5 (2) 

6 . Suz 
2.S6(3) 
Sz (8) .3 
3. J3 

3. U6 (2) 
2F4 (2) 

2 . Ru 
S8 (3) 
M2 3 
3. ON 
J4 



Irred. character no. 



10 

2 
29 

2 
49 
19 
25 
22 

3 

153 
2 
4 

22 

47 
2 

37 
2 
3 

31 
2 



TABLE II: Generators of a subgroup G of Sp(3, 2) of order 2(d 2 - 1) = 160. T he group G acts transitively on the non-zero elements of the 
phase space V — F| and induces a unitary design, as described in Section llV C| 



( 2 2 2 \ 

12 2 

12 2 

\ 1 1 / 



/ 2 1 ^ 

2 

10 2 

\ 2 / 



/ 2 ^ 

2 

2 2 

\ 1 2 / 



/ 1 1 ^ 

10 2 

1 

\ 1 / 



[ 1 0^ 

2 

2 

V 1 / 



^ 2 ^ 

2 

2 

V 2 / 



TABLE III: Generators of a matrix group G of order 120. Up to a similarity transformation, the group gives rise to a unitary 5-design with 60 
elements. In fact, G is an irreducible representation of SL(2, Fs) affording the character X.3 as listed in the GAP character library |33|| . The 
repersentation has been constructed using the "Repsn" package 13411 . Below, uj = exp(^-i). 

1 \ f -LU 11 - LU 14 -LU 11 - LU 14 \ f -LU - LU 2 - U 4 - LU S - 2{LU 
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